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ABSTRACT 

Aims. We study the two main constituent galaxies of a constrained simulation of the Local Group as candidates for the Milky Way 
(MW) and Andromeda (M31). We focus on the formation of the stellar discs and its relation to the formation of the group as a rich 
system with two massive galaxies, and investigate the effects of mergers and accretion as drivers of morphological transformations. 
We also assess the effects of varying the assumed feedback model on our results by running two different simulations, a first one 
where only supernova feedback is included and a second where we additionally model radiation pressure from stars. 

Methods. We use a state-of-the-art hydrodynamical code which includes star formation, feedback and chemical enrichment to carry 
out our study. We use our two simulations, where we include or neglect the effects of radiation pressure from stars, to investigate the 
impact of this process on the morphologies and star formation rates of the simulated galaxies. 

Results. We find that the simulated M31 and MW have different formation histories, even though both inhabit, at z = 0, the same 
environment. These differences directly translate into and explain variations in their star formation rates, in-situ fractions and final 
morphologies. The simulated M31 candidate has an active merger history, as a result of which its stellar disc is unable to survive 
unaffected until the present time. In contrast, the MW candidate has a smoother history with no major mergers at late times, and forms 
a disc that grows steadily; at z = 0 the simulated MW has an extended, rotationally-supported disc which is dominant over the bulge. 
Our two feedback implementations predict similar evolution of the galaxies and their discs, although some variations are detected, the 
most important of which is the formation time of the discs: in the model with weaker/stronger feedback the discs form earlier/later. In 
summary, by comparing the formation histories of the two simulated galaxies, we conclude that the particular merger/accretion history 
of a galaxy rather than its environment at the LG-scales is the main driver of the formation and subsequent growth or destruction of 
galaxy discs. 
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1. Introduction 


The formation of late-type spiral galaxies in the context of the 
A-Cold Dark Matter (ACDM) cosmological framework is a long 
standing problem in astrophysics. In hierarchical models such as 
ACDM, the presence or absence of a strong, rotationally sup¬ 
ported stellar component is influenced by processes acting at 
many scales as galaxies evolve; baryons condense in the centres 
of dark matter haloes whilst retaining a considerabl e fraction of 
their angular momentum set at large scales (Fall & Efst athioul 
Il980h and a combination of viscous and dissipative processes 
in the gaseous component allows the formation of a disk from 
which a fraction is converted into a stellar component. Thicker 
stellar comp onents result from dynamica l processes such as disc 
instabilities (I Sell wood & Binnev, 2002) and the close encoun¬ 
ters and mergers (e.g. lOuinn. Hernquist & Fullagaii 1 19931) that 
are ubiquitous in hierarchical growth mode ls suc h as ACDM, 
though some simulations (e.g. [Bird et all 12013 ) alternatively 
suggest the older stellar disks may have been born thick. 

Numerical simulations are the best tool to study galaxy for¬ 
mation within ACDM, and to investigate the relationship be¬ 
tween the formation of stellar discs and the hierarchica l build¬ 
up of the haloes containing them (iScannapieco et all [2009; 
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ISales et all 120121) . During recent years, important progress has 
been made in the field, and simulations currently appear to be 
reaching the level of sophistication required to study the de¬ 
tails of the assembly of galaxies and their stellar discs. In fact, 
while a first generation of models inevitably suffered from the 
angular momentum catastr ophe (e.g., iNavarro & Benzl 1 19911: 
ISteinmetz & Navairo[ll999l) . more recent codes showed that var¬ 
ious types of feedback provide galaxies with self-regulated star 
formation activity, in turn allowing the formation of young , 
rotationally-supported stellar discs (IScannapieco et all [20121) . 
Feedback processes produce a decoupling between the forma¬ 
tion of a gas disc after its condensation within a halo and the 
onset of star formation activity; in this way, stellar discs can 
form later, in more quiescent periods rather than during the 
early, violent times characteristic of the formation of the haloes. 
When discs form later, they also have a greater chance to sur¬ 
vive a nd can grow to form exten ded, highly-rotating systems by 
z = 0 (IScannapieco et al.. 2009). At the heart of the problem of 
disc formation/survival is therefore the problem of understand¬ 
ing how star formation proceeds in galaxies, and what are the 
possible regulating mechanisms and their relative efficiencies. 

For Milky Way-mass galaxies, the most efficient feedback 
mechanism is believed to be that produced by stars, via the 
injection of energy into the interstellar medium during differ¬ 
ent phases of stellar evolution. In particular, various models for 
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feedback from supernovae (SNe) have been included in simu¬ 
lations in the last decade, and demonstrated that SN feedback 
is indeed capable of regulating star formation and of produ cing 
stella r discs with high angular momenta (e.g.. lOkamoto et al~ . 
200 51: iGovernato et~aTI. 120071: IScannapieco et all I2008L 12009 : 
Piontek & Steinmetz, 201 1]). However, stellar discs in these 


models are usually subdominant over the massive bulges that 
form early. The more recent inclusion of additional feedback 
in the form of radiation pressure produced by massive stars 
showed that this process provides significant kinetic energy to 
the gas, particularly at early times, helping to produce even 
younger, more massive galaxy discs sim i lar to those observed 
(lAumer et all 120131: IStinson et all 120131: iHopkins et all 12014 
lAgertz & Kravtsovll2014l). 

Des pite the recent progress, the Aquila code comparison 
project dScannapieco et all 12012b showed that, even when the 
formation history of a halo is fully specified, different models 
predict the formation of a galaxy with different gas fractions, 
morphologies, stellar masses and angular momentum content. 
Aquila clearly demonstrated that the primary uncertainty in these 
types of simulations is related to the processes of the unresolved 
baryonic physics that result in star formation and feedback and 
their respective modeling, and generated consensus on the need 
to produce similar results using different codes, despite their 
individual successes. Efforts to overcome these problems are in¬ 
deed ongoing, both in terms of direct comparison between codes 
dKim et all 12014) and of the development of more physically- 
motivated routine s to tre at star formati o n and fe edback (e.g.. 


Christensen et al.l 20121: Agertz et al. 

20 lH lAumer et al. 

2013 

; Creasev, Theuns & Bower 

20lit Stinson et al. 

2013 

: Vogelsberger et al. 2013: Agertz & Kravtsov 2014 

Hopkins et al. 2014 Muranteetal. 2015t Truiillo-Gomez et al.1 


20151) . 


An important motivation to study the formation of late-type 
spiral galaxies is the membership of our own galaxy, the Milky 
Way, within that family. The MW provides an ideal place to con¬ 
strain galaxy formation theories, as we have detailed observa¬ 
tional data on a large number of stars residing in its different 
stellar components. A salient characteristic of the MW is that it 
lives in a group environment, known as the “Local Group” (LG), 
with another large spiral, Andromeda (M31), lying at less than a 
Mpc distance from it. A number of observational studies showed 
that the immediate environment of galaxies indeed plays a role in 
the determination of their p roperties at the s cales of galaxy clus¬ 
ters and rich groups (e.g., Dressier, 1980; Herm it et all 119961: 
iGuzzo et aLlll997l:lBlanton et all 120031:iGirardi et alll2003l) . al¬ 

though for systems like the LG this i s more controversial (e.g., 
iBaheet all 120131: IZiparo et all l2013h . On the other hand, sim¬ 
ulation studies have almost exclusively used initial conditions 
(ICs) that follow the formation of a galaxy which is isolated from 
massive neighbours at the Mpc scales (see, however. iFew et all 
120121) . A study on possible effects of environment on the forma¬ 
tion of the MW and M31 is therefore required in order to under¬ 
stand if these are representative of galaxies of a similar mass and 
to better interpret observational results. 

In this work, we use constrained ICs of the Local Universe, 
where a LG-like pair of galaxies form. We study the formation 
and evolution of these two LG-galaxies, candidates for the MW 
and M31, using a state-of-the-art hydrodynamical code for cos¬ 
mological simulations. We use these simulations to investigate 
how the LG galaxies could have formed, and to look for possi¬ 
ble imprints of the formation of the LG on their final properties. 
In particular, we focus on the formation of the main stellar com¬ 
ponents of the galaxies, and investigate possible correlations be¬ 


tween the formation of stellar discs and large-scale effects such 
as merger events and misalignment between the gaseous and 
stellar discs. The main questions we want to address here are 
whether our simulation predicts the formation of two large spi¬ 
ral galaxies similar to the MW and M31 and what are the mech¬ 
anisms determining/affecting their formation. We will show that 
only the simulated MW has a stable stellar disc that survives to 
z = 0; while in M31 various mergers continuously affect the 
disc, and it produces a galaxy with a much lower level of rota¬ 
tional support. We find similar results using two different ver¬ 
sions of the code that we use to test whether results are affected 
by the particular treatment of feedback pr ocesses. It is wo rth 
notin g that, as shown in the Aquila Project (IScannapieco et all 
120121) . the final morphology of a galaxy depends sensitively on 
the amount of star formation during the early epochs of the for¬ 
mation of galaxies, and therefore on the treatment of feedback 
processes which regulate the star formation activity. 

Other studies using our simulations are presented in a series 
of companion papers. iNuza et al.1 ( 2014bl) analysed the proper¬ 
ties of the neutral, warm and hot gas within and around the 
haloes of the MW and M31 candidates in one of our simulations, 
including accretion/ejection rates of the different components 
and HI covering fractio ns that we contrasted to observational re¬ 
sults. In lCreasev et al.1 (2015 ). we further investigate the relation 
between environment and galaxy properties, by comparing the 
properties of our simulated MW and M31 with those of galaxies 
of similar mass but formed in more isolated environments. Our 
findings suggest that some of the properties of galaxies, such as 
the star formation rates, may indeed be directly affected by the 
environment at the Mpc scales; while others, like the disc stellar 
mass fraction, are more sensitive to the merger and approach of 
satellite systems to the inner regions of the galaxies. (For other 
simulation studies analysing the relation between galaxy proper¬ 
ties and environment in constrained simulations - albeit at much 
larger scales - see e.g. INuza. Dolag & Sarol 120101 : INuza et all 
120143) . 

The outline of our paper is as follows: in Section |2] we 
describe the initial conditions and code used for our study; 
Section [3] describes the assembly of the LG and the LG galax¬ 
ies in our simulations; and in Section 0] we discuss the morpho¬ 
logical transformations of the simulated galaxies during their 
evolution, and investigate the relation between the evolution of 
morphologies and the merger histories. Finally, in Section 0 we 
summarise our results. 


2. The simulations 

We use for this study a simulation designed to produce a LG-like 
system in the right cosmological environment, and investigate 
the formation of its two main constituent galaxies, candidates 
for the Milky Way and Andromeda. In this section we describe 
the numerical methods needed to carry out such a study, and 
summarize the numerical set-up of the simulations. 

2.1. Simulation code 

We use the Tre e-PM Smoothed Particle Hydrodynam ics (SPH) 
code gadget3 dSpringeli 120051: ISpringel et all 120081) . and run 
two different simulations where additional updates are in¬ 
cluded. The first simulation, labelled “CS” throughou t this 
paper, uses the extensio ns of IScannapieco et al.l (120051) and 
IScannapieco et al.l (120061) for star formation, metal-dependent 
cooling, chemical enrichment and (thermal) feedback from 
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Type II and Type la supernova (SN) explosions, a multi¬ 
phase model_Joi_the_gas__component, and a UV background 
field (lHaardt & MadaulT 19961) . The second simulation, labelled 
“MA”, additionally includes the effects of radiation pressure 
from stars and kinetic feed back from SN explosions, as de¬ 
scribed in lAumer et al.1 (1201 3). 

We have shown in previous work that the CS model is able 
to reproduce the formation of galaxy discs from cosmological 
initial conditions, alleviating the angular momentum problem. 
However, although discs have realistic sizes and angular mo¬ 
mentum content, overly massive bulges are also fo rmed which 
are usually dominant (in mass) over the discs (IScannapieco et al.l 
1200 8ll2009i 120 1 0ll20 1 lL 120 1 2l) . The more recent MA implemen¬ 
tation is an update to the CS model. In this case, early stellar 
feedback from new born stars and kinetic feedback from SN 
explosions is additionally considered, which leads to stronger 
effects particularly at high redshifts. The MA model produces 
galaxies with a high angular momentum content, moderate stel¬ 
lar masses and dominant discs, in better agreement with obser¬ 
vations of spiral galaxies. It is therefore relevant for our work to 
compare the results obtained with the two models, which allows 
us to test the robustness of our results, particularly in relation to 
the effects of mergers during the evolution of galaxies. 

In both models, star formation takes place in dense (Ah > 
0.03 cm -3 ) regions of converging flow, with a star formation rate 
per unit volume equal to 

P* = c, (1) 

Ulyn 

where p and Td yn = 1 / ^4nGp are respectively the density and 
dynamical time of the gas particle, and c* a star formation effi¬ 
ciency (that we set to c * = 0.1). Our simulations u se the mul¬ 
tiphase model described in IScannapieco et all (120061) . where the 
neighbour list of any particle i with entropy A t is constructed ig¬ 
noring neighbours j with Aj < 50 A t (provided they are not part 
of a shock to avoid unphysical behaviour). In this way, a multi¬ 
phase medium is generated, where particles of varying densities 
and temperatures (i.e. cold and hot ) can coexist in the same spa¬ 
tial region. These cold and hot phases are however not fixed, as 
the selection of particles to be ignored as neighbours depends 
on the relative properties of each pair of particles i, j and is not 
ba sed on any pre-fixed value s for different phases. As explained 
in IScannapieco et al . (12006 ), not only this model allows coex¬ 
istence between a dense and a diffuse phase, but also a more 
efficient deposition of the supernova energy. 

In the CS model, each star particle can “explode” a maxi¬ 
mum of two times, first as SN Type II (SNII) and later on as SN 
Type la (SNIa). During each SN explosion, the star distributes 
chemical elements (according to the corresponding SNII and 
SNIa chemical yields) and energy into neighbouring gas parti¬ 
cles. The energy and chemical production of each exploding star 
is distributed in equal proportions to its hot (non star-forming) 
and cold (star-forming) gas neighbours. We assume a canonical 
0.7 x 10 51 ergs of energy to be released per SN. The distribu¬ 
tion of chemical elements occurs at the time of the explosion, 
as well as the release of energy to the hot phase. However, for 
the cold phase, feedback energy is accumulated in a reservoir, 
until it is sufficient to allow the particle thermalize with the lo¬ 
cal hot phase. In practice, this works as follows: for each cold 
gas particle, we calculate the average density and entropy over 
its hot neighbours (those particles that were ignored in its neigh¬ 
bour list according to our multiphase model) and calculate the 
energy needed for the cold particle to reach these values, A E. 


This energy is then compared to the energy stored in the reser¬ 
voir E ygs , and when E res > A E the cold particle gets the reservoir 
energy as thermal energy. Note that the value of A E is different 
and time-dependent for each cold particle, which makes the code 
suitable to adapt to the cold/hot regions with different tempera¬ 
tures and densities, which are present at different times during 
the formation of a galaxy and at different spatial regions. 

The MA model works similarly to the CS model; however, 
in this case, the SN energy is split into a kinetic and a thermal 
component. The kinetic energy transferred from the exploding 
star to the gas neighbours is calculated assuming that the mo¬ 
mentum of the SN ejecta is conserved, and characterized by an 
outflow velocity L> out : A p = Am • y ou t. We use y out = 3000 km s -1 , 
with momentum given in the radial direction. As the momen¬ 
tum is distributed into neighbouring particles, which typically 
translates into a velocity kick per particle of Av{ ~ 20- 30 km 
s -1 . The remaining energy (i.e. the total energy released per ex¬ 
ploding star minus the kinetic energy that has been already been 
transferred to the neighbouring gas particles) is stored into the 
thermal reservoir of the gas receiving feedback, and is treated 
in the same way as in the CS model. The MA model includes 
a treatment for radiation pressure, where the rate of momentum 
deposition in the gas is parametrized as 


Pip = (1 +tir) - (2) 

c 

with L(t) being the UV-luminosity of the stellar population and 
tir an infrared optical depth which depends on the density, 
metallicity and velocity dispersion of the gas around the star. 
This last dependency makes the effects of radiation pressure 
larger in systems with higher velocity dispersions, typical of 
high-z galaxies, allowing a more efficient regulation of star for¬ 
mation at early epochs and higher star formation levels at late 
times. 

Note that, as described in lAumer et al.l ( 20131 ) . the MA model 
considers additional updates (new treatment of chemical enrich¬ 
ment and metal diffusion, AGB stars, and element-by-element 
cooling) that we do not include here, and therefore the only dif¬ 
ferences between our two simulations is the inclusion of radi¬ 
ation pressure and the additional kinetic feedback assumed for 
SNe. In order to better compare results from our two simula¬ 
tions, we have also used the exact same input parameters for star 
formation and feedback when appropriate (which in some cases 
differ from those assumed in lAumer et al . 1120131 ) . and also used 
the same choice for the gravitational softening, which was (for 
the high resolution particles) 0.7 kpc, fixed in physical coordi¬ 
nates since z = 3 and fixed in comoving coordinates at earlier 
times. In this way, any difference in the results can be directly 
associated to the effects due to the radiation pressure of young 
stars and kinetic feedback during the SN phase. 

We note that some properties of the simulated galaxies, par¬ 
ticularly the detaile d final morphology, are e xpected to be sen¬ 
sitive to resolution (IScannapie co et al., 2011) and to the choice 
of the gravitational softening (lAumer et al . 1120131 ) . However, al¬ 
though higher resolution or/and lower values for the gravita¬ 
tional softening will allow the vertical/radial structure of discs 
to be better resolved, the presence or absence of a disc will not 
be affected by these choices (as long as a minimum resolution is 
considered). Our current work focuses on the evolution of discs 
and therefore our results are still valid, even when the spatial res¬ 
olution is not high enough to resolve the radial/vertical structure 
of discs with a large number of resolution elements. 
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2.2. Initial conditions and simulation set-up 

The initial conditions (ICs) used for this study are con¬ 
strained initial conditions of the Local Universe, that are part 
of the CLUES (Constrained Local UniversE Simulation^]) 
project. The ICs reproduce, by construction, the known 
dynamical properties_ of our loc al environment at the 
present-day (see GotUober, Hoffman & Yepesl 120101 and 
lYepes. Gottlober & Hoffmanl 120141 for details) and use the 
zoom-in technique, in a periodic box of 64 h ~ l Mpc on a side. 
In the centre of this cube, a high-resolution region of 2 h ~ l Mpc 
comoving radius contains dark matter and gas particles, with a 
mass resolution of 2.8 x 10 4 5 6 M© and 5.6 x 10 5 M© respectively. 
Outside the high-resolution region, the mass distribution is 
described by low-resolution collisionless particles, as usual in 
zoomed simulations. 

The ICs are consistent with a ACDM universe with WMAP- 
5 cosmological parameters: Qm = 0.279 (matter density), Qa = 
0.721 (dark energy density), Qbar = 0.046 (baryon density), 
Ho = 100 h km s -1 Mpc -1 with h = 0.7 (Hubble parameter), and 
o ~8 =0.8 (normalization of the power spectrum), and correspond 
to a starting redshift of z = 50. 


2.3. Milky Way and Adromeda candidates 


At z = 0, the simulated LG has two main constituent galaxies, 
with total masses similar to the MW, that we take as candidates 
for the MW and AndromedcB We will refer to these galaxies 
as G1 (candidate for M31) and G2 (candidate for MW0 and, 
where appropriate, we will add an acronym CS/MA to denote 
the simulation run. The assignment of our galaxy candidates to 
the real counterparts is based on their mass, G1 being the most 
massive at z = 0 (see Table [T]). Additionally, with this choice, 
we recover a similar relative ori enta tion between the galaxies in 
comparison to observations dNuza et all 12014bh . 

Table \T\ shows the main z = 0 properties of G1 and G2 in our 
two runs: R 200 , defined as the radius where the density contrast 
is 200 times the critical density, and the total (M 200 ), gaseous 
(Afgas) and stellar (M stars ) masses, all within /?200- The present- 
day masses of G1 (~1.7 xlO 12 M 0 )andG2 (~1.2xl0 12 M©) are 
similar to the values estimated for Andromeda (- 1.4 x 10 12 M o ; 


Way (~10 12 M«; Will 

kinson & Evans 1999j; 

Smith et al. 

2007; 

Watkins. Evans & An 

201(llBow et al. 2012 

Kafle et alj 

2012; 

Piffl et al. 2014) from observational studies. 


3. The assembly of the LG and the LG galaxies 

The assembly history of the LG and of the LG galaxies (Gl, 
candidate for M31 and G2, candidate for MW) is seen in Fig. [U 
where we show the dark matter, gaseous and stellar mass dis¬ 
tributions (in an arbitrary projection) at z = 2,1,0.5 and 0. The 
different panels of Fig. [T] are centred at the centre of mass of the 
G1-G2 system at the corresponding cosmic epochs, and the two 
galaxies are surrounded by solid (Gl) and dashed (G2) circles 
indicating their R 200 radii. At each redshift, the plots show the 


Table 1. Main properties of the Gl and G2 galaxies in our two 
simulations at z = 0: radius R 200 , and total, gaseous and stellar 
(M 200 , M gas and M stars ) masses within within T^oo- 


Galaxy- 

Model 

R 200 

[kpc] 

M 200 

[10 10 Mo] 

^gas 

[IO 10 Me] 

^Vtars 

[10 10 Me] 

Gl-CS 

244.9 

167.9 

6.67 

7.95 

Gl-MA 

239.1 

155.1 

6.13 

6.76 

G2-CS 

222.2 

125.1 

5.69 

6.24 

G2-MA 

219.5 

120.6 

6.39 

5.48 


mass distributions in a cubic box of 2.5 cMpcQ side length that 
corresponds to the high-resolution region. The different colours 
span 4 orders of magnitude in logarithmic scale, and we adopted 
the same colour scale at all redshifts in order to highlight time 
variations in the distribution^. 

At z = 2, the characteristic filamentary structure at Mpc 
scales is clearly seen; Gl and G2 are 550 kpc apart (i.e. more 
than a Mpc apart in comoving coordinates), and belong to differ¬ 
ent filaments. The Local Group evolves quickly and, by z = 1, 
defines a system where the mass distributions of its two main 
galaxies are interconnected - a feature that is particularly visi¬ 
ble in the gas component - and have collapsed to the same fila¬ 
ment. The two galaxies continue separating from each other until 
z - 0.5 when they reach a maximum separation of - 900 kpc. 
After this time, Gl and G2 start approaching each other, to reach 
a final separation of about 800 kpc. The stars are located within 
the very centres of the dark matter haloes, resulting in much less 
extended distributions compared to the gaseous and dark matter 
counterparts. A significant number of smaller galaxies and satel¬ 
lites of Gl and G2 are observed, particularly at high redshift; 
most of them have been accreted by the Gl and G2 haloes by 
z = 0. 

Although Gl and G2 are, at z = 0, both part of the simu¬ 
lated LG, their respective growth history are somewhat differ¬ 
ent, as evidenced from the evolution of their total masses and 
radii shown in the left-hand panels of Fig. [2] In the case of Gl, 
the size and mass increase rapidly from z = 3 to z = 0, particu¬ 
larly after z ~ 1 owing to the accretion of a number of satellite 
galaxie^. These satellites are being accreted at approximately 
the same time around z = 0.9 - 0.7, and contribute a significant 
amount of mass to Gl, determining the increase and variations 
of /?200- Slightly after z = 0.5, the satellite galaxies are in their 
last apocentre, and the evolution of Gl becomes more quiet until 
the present time. In contrast, G2 has a more steady evolution at 
all times, with the most significant changes at z - 1-1.5, which 
also result from the entrance of a relatively small satellite galaxy. 
The agreement between the evolution of Gl and G2 in our two 
simulations is excellent in terms of the masses and radii. 

The different growth of the haloes of Gl and G2 also re¬ 
flect in the assembly of the baryonic components, with a much 
smoother evolution of G2 compared to Gl, as shown in the 
right-hand panels of Fig. [2] However, as the evolution of the 


1 http://www.clues-project.org/ 

2 Note that, as the constraints can only be imposed on large scales, 
the properties of the two main galaxies will not necessarily resemble 
those of the actual MW and M31 in detail. 

3 Note that in iNuza etal. ( 2014b h we instead refer to Gl and G2 as 
M31 c and MW C respectively. 


4 For short, throughout the paper we will use a prefix c to length units 
to denote the use of comoving coordinates. 

5 These figures correspond to our CS run although, at these scales, 
the results are qualitatively the same as those obtained with the MA 
model. 

6 We discuss the merger histories of Gl and G2 in Section[4] 
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Fig. 1. Spatial distribution (in an arbitrary projection) of the dark matter (upper panels), gas (middle panels) and stars (lower panels) 
in our simulation, in a cube of 2.5 cMpc on a side and for redshifts z — 2,1,0.5 and 0. The R 200 of G1 (solid line) and G2 (dashed 
line) at each redshift are indicated, and the colour scale (of each mass component) is fixed in order to highlight differences in the 
distributions at different times. The plots correspond to run CS, but results at these scales are very similar for both simulations. 


gaseous/stellar masses are affected not only by the growth of the 
haloes with time, but also by the transformation of gas into stars 
and the return of mass from the stars to the gas owing to stellar 
evolution, larger differences between the two runs are found for 
the baryonic component^ As a result of the stronger feedback 
assumed in model MA, galaxies have in general lower stellar 
and higher gaseous masses in this run, compared to simulation 
CS. This effect is already noticeable at high redshift, due to the 
early feedback associated to massive, short-lived stars, that acts 
shortly after the onset of star formation activity in the galaxies. 

The differences in stellar/gas masses between results of runs 
CS and MA are a consequence of the different star formation 
rates (SFRs), shown in the upper panels of Fig. [3] To calculate 
the SFRs, we considered all stars that end up in the main galaxy 
at z = 0 (i.e. within 30 kpc), and therefore include stars formed 
in-situ (within the main progenitor at each time) and in satellite 
systems that were later accreted. The SFRs are higher at early 
times and gently decline at more recent epochs. Feedback effects 
start before in run MA compared to run CS, as radiation pressure 


7 Note, however, that the simulations presented here only consider 
mass return from SN explosions, while they do not include a treatment 
for the mass return of intermediate-mass stars that would result in a 
larger correction. 


acts right after the onset of star formation. In fact, the first star 
formation burst produced in run CS (in both galaxies) is absent 
in simulation MA. In the case of Gl, we find that the SFRs in 
model MA are in general lower compared to model CS, with 
the most important differences at very early times and around 
6-8 Gyr. For G2, differences are significant (additionally to the 
very early epochs) at 4 - 6 Gyr, with lower SFRs for model MA, 
and after 9 Gyr of evolution, with higher SFRs for model MA 
compared to CS. 

The lower panels of Fig. [3] show the corresponding in-situ 
fractions, /in-situ, which help us to understand variations between 
the evolution of Gl and G2 that result from the different SFRs 
of the progenitor galaxies and their satellites. An important dif¬ 
ference between Gl and G2 (in both runs) shows up very clearly 
in these plots: Gl has much lower in-situ fractions at all times 
compared to G2. Only at low redshifts does Gl have /in-situ ~ 1, 
while at earlier times there is a high contribution of stars that 
form in other systems than the main progenitor. In contrast, G2 
has higher in-situ fractions at all times: /in-situ ~ 1 for at least 
the last 6-8 Gyr. As we discuss in Section 0 these differences 
can explain the origin of the final morphologies of the simulated 
galaxies, as Gl has a much more active merger history which 
makes its stellar disc much more susceptible to partial or total 
destruction. 
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Fig. 2. The size/mass evolution of G1 (solid line) and G2 (dashed line) as a function of redshift. We show the evolution of R 200 
(upper-left panels), the total mass (M 200 , lower-left panels), the gas mass (M gas , upper-right panels) and the stellar mass (M star , 
lower-right panels), all within T^oo- For Gl, the peak seen in all panels at z ~ 0.65 is owing to the emergence of two large objects 
within R 200 (see further discussion in Section 0]). 



Fig. 3. Star formation rate of Gl (left-hand panel) and G2 (right-hand panel) as a function of cosmic time since the Big Bang for 
simulation runs CS (solid lines) and MA (dotted lines). The lower panels show the corresponding in-situ fractions, i.e. the stellar 
mass fraction formed in the main progenitor at each time interval. 


4. Galaxy morphologies and their relation to the 
formation/merger history 

The figures of the previous section already revealed that Gl and 
G2 exhibit differences in their growth/star formation history. In 
this section, we investigate how these differences are translated 


into variations in the final properties of the galaxies, particularly 
in their morphologies/dynamical state. 

The stellar distributions of the central, main component of 
Gl and G2 in our two runs and at z = 0 are shown in Fig. 0] 
for face-on and edge-on views, in an area of (40 kpc) 2 . To cre¬ 
ate these maps, we have rotated the galaxies such that the stel- 
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Fig. 4. Projected stellar mass density of G1 (left-hand panels) and G2 (right-hand panels) showing edge-on and face-on views for 
simulation runs CS (upper panels) and MA (lower panels). Also shown are the velocity field (arrows) and isodensity contours (solid 
lines) in the corresponding plane. 




Fig. 5. Distribution of stellar circularities of G1 (left-hand panel) and G2 (right-hand panel) for simulation runs CS (solid lines) and 
MA (dotted lines) at z = 0. 


lar angular momentum within a radius of 30 ckpc lies in the 
z-direction. We also plot isodensity contours at fixed density 
values in order to highlight differences in the mass distributions 
of the two simulations. 

Immediately apparent is that, in both runs, G1 is a more 
spheroidal galaxy than G2, which has a well-defined, extended, 
disc-like component. The presence/absence of a disc can be in¬ 
ferred from the projected velocity fields overplotted to the stel¬ 
lar densities. These show that G1 has a certain degree of rota¬ 
tional support that is however not very significant. Compared to 
Gl-CS, the rotation in Gl-MA is much more coherent, as evi¬ 
denced not only by the tangential velocities in the face-on view 
but also by the low velocities observed in the edge-on projection. 
A much higher degree of rotation is found for G2 in both simula¬ 
tions, where a well-defined, extended disc component is present. 


As shown later (Figs. [5] and [7]) neither G2-CS nor G2-MA has 
a dominant, massive bulge in the centre, and neither shows evi¬ 
dence of a bar. Note, however, that due to its younger age (Fig. 0 
see below), the disc of G2-MA is thinner compared to that of G2- 
CS. Also note that, both for G1 and G2, we find that the stellar 
haloes are more extended in run CS compared to run MA, which 
is a direct consequence of the stronger feedback assumed in the 
latter which prevents the formation of the old stars that make up 
this component. 

In order to better quantify the dynamical state of G1 and G2 
and compare results from our two simulations, we show in Fig. [5] 
the corresponding distributions of stellar circularit i es. Th e circu¬ 
larity of a star i is defined as in lScannapieco et all (l2QQ9h : 


7 





























































Scannapieco et al.: A constrained simulation of the Local Universe 


where j\ is the angular momentum of the i star particle perpen¬ 
dicular to the disc plane (i.e. in the z— direction), and j c ^ c the 
angular momentum expected for a circular orbit at the star’s ra¬ 
dius. Circularities of the order of 1 are expected for rotationally- 
supported particles while spheroidal components form more or 
less symmetric distributions around 6 = 0 (non-rotating bulges), 
0 < 6 < 0.5 (rotating bulges) or -0.5 < e < 0 (counter-rotating 
bulges). 

As inferred from Fig.0 Fig. [5] confirms that G1 has an im¬ 
portant spheroidal-like component. In Gl-CS, the spheroid has 
an overall rotation (the peak of the distribution is at positive e), 
while in Gl-MA the spheroid does not rotate. Furthermore, Gl- 
MA has a disc-like component in rotational support. As we show 
later, the presence/absence of the disc in Gl-CS relates to the 
fact that this galaxy has experienced, before z = 0, a number of 
mergers which destroyed a pre-existing disc, while these events 
have not yet happened in Gl-MA at z = 0. In the case of G2, the 
presence of a well-defined, rotationally-supported disc is clear in 
both runs. In addition, G2-CS has a non-rotating, distinguishable 
although not-dominant bulge. 

The SFRs of G1 and G2 in our two simulations, shown in the 
previous section (Fig. [3]) and, in particular, the corresponding in- 
situ fractions, already gave some indications of what the final 
morphology of the galaxy would be. In Gl, star formation is ex¬ 
tremely bursty, and the in-situ fractions are relatively low, which 
is indicative of a large contribution of stars from satellite galax¬ 
ies. In general such a contribution from satellites will reduce the 
disc fraction, unless the angular momentum of the satellite is 
very closely aligned with that of the progenitor galaxy. In the 
case of G2, the smoother SFRs together with very high in-situ 
fractions, particularly at low redshifts, suggest the formation of 
a galaxy with a much more significant disc which did not suffer 
disruption from violent events. 

In the rest of this section we show that, in fact, Gl has a much 
more violent merger history than G2 which has a much smoother 
evolution. It is worth noting that the in-situ fractions of our sim¬ 
ulated galaxies are similar in runs CS and MA which is an indi¬ 
cation that our results on the final morphologies (particularly of 
Gl) are robust and not triggered by the assumed feedback model. 
Mergers seem to be as disruptive in both runs, regardless differ¬ 
ences in the exact merging times of the satellites (see below) and 
the different final stellar masses of the galaxies. 

Fig-0 shows the merger histories of Gl and G2 in our two 
simulations, which we constructed by iden tification o f haloes us¬ 
ing the substructure finder code subfind (ISpringel et all 1200ll) 
and then matched between snapshots with a metric that weights 
intersecting particles by the ranking of their binding energy. We 
identified all mergers with a mass ratio larger than 0.1 % in terms 
of their total mass, and denoted any merger greater than 30% as 
a ‘major’ merger and any greater than 5% as an ‘intermediate’ 
merger. Table [2] describes statistics - redshift and mass ratio - of 
the intermediate and major mergers as the satellites fall within 
the R 200 of the progenitor galaxy. We also show the mass ratios 
in terms of the stellar mass, as these might be more indicative of 
the importance of mergers in the central regions. 

Whilst the accretion of a progenitor into R 200 will clearly 
play a role in the evolution of the galaxy, the events at much 
smaller radii have a more direct impact, as an infalling object 
can orbit the galaxy for several Gyrs (or more) before merging 
into the central object. For this reason we followed the merg¬ 
ing subhaloes also after their entrance into R 200 , and calculated 
the times and mass fractions when they enter the inner 30 ckpc, 
where we study the stellar and gaseous components. This ra¬ 
dius is a good choice as it is both large enough to encompass 


the galaxy (see Fig. n but also small enough that objects whose 
centres have passed this close will significantly interact. For both 
simulations, the last three columns of Table [2] show the redshift 
at which satellites cross the 30 ckpc threshold, as well as the stel¬ 
lar mass ratio at that time. Note that, in order to match results of 
runs CS and MA, we followed the mergers of satellite galaxies 
even though, in some cases, they have not yet merged/disrupted 
with the main progenitor and are still identifiable, independent 
systems at the present time. These are indicated in the final col¬ 
umn of the Tabled with annotation ‘sat’. Also note that the times 
of the mergers, particularly when the satellites reach the inner 30 
kpc, are different in the two runs, as well as the merger ratios. 

From Fig. [6] and Table [2| we can see that, as inferred from 
the growth of the haloes and the in-situ fractions, Gl has a more 
active history than G2, with 6 major/intermediate mass mergers, 
compared to only 4 in the case of G2. Perhaps more significant 
is that at late times (z < 1) Gl accretes three very large objects 
indeed, while for G2 there are no mergers in this period (consid¬ 
ering their entrance to R 200 X and its only mergers for z e [0, 3] 
are at much more modest stellar fractions. Although similar in 
general terms, the details of the merging process of Gl and G2 
in our two runfl show some important differences. In particular, 
in Gl-CS only one (and small) satellite survives at z = 0 as an 
independent object, while the same object has been already dis¬ 
rupted in run MA. More important is the fact that the two satel¬ 
lites which merge later (denoted Gle/f) have been disrupted by 
z = 0 in run CS but not yet in run MA. For G2, we find less varia¬ 
tion between the simulations and, in any case, the merger events 
of G2 are at much lower fractions and occur earlier compared 
to those of Gl. For this reason, it is expected that the merger 
history will have a much stronger impact in the formation of Gl 
compared to G2. 

The growth/destruction of galaxy discs is intimately related 
to the merger histories of galaxies. In fact, we find a clear rela¬ 
tion between periods of disc growth/destruction and the merger 
history of our simulated galaxies. Fig. [7] shows the evolution of 
/disc, a measure of the disc-to-total (D/T) ratio, for Gl and G2 
in runs CS and MA. Note that, from Fig. 0 it is clear that the 
distribution of circularities is a very clear indicator of the pres¬ 
ence or absence of a disc-like stellar component and thus a single 
quantitative measure of this is the fraction of stellar mass with 
circularity largei0 than 0.5, 

/disc =f(e> 0.5). (4) 


Note, however, that /di SC can be non zero even if there is no 
disc, if the distribution of circula rities is sufficien tly broad. On 
the other hand, previous work ( Scannapi eco et ap2010l: see also 
lAbadi et aDl2003L iGovernato et al.ll2007l) showed that kinematic 
D/T decompositions, as the one we present here, give ratios 
that are systematically and significantly lower than those ob¬ 
tained with photometric decompositi ons (i.e. comparab l e with 
observed D/T ratios). For example, in lScannapieco et all ( 2010 ) 
we showed that D/T can increase from 0.2 (in a kinematic de¬ 
composition) up to 0.5-0.7 (in a photometric decomposition). In 
any case, as the primary focus of our analysis is to find an appro¬ 
priate measure of the evolution of a galaxy’s morphology (and 
not to make a detailed comparison of simulated and observed 
D/T ratios), we can safely use /di SC as it sufficiently captures the 
prominence of the disc as a function of time in our galaxies. 


8 Note that this can explain some of the differences in the properties 
of Gl and G2 in runs CS and MA, e.g. morphologies, that we discussed 
in this and the previous sections. 

9 In the calculation of / disc , we consider only particles of the main 
subhalo, excluding satellites. 
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Fig. 6. The merger history of G1 (left-hand panels) and G2 (right-hand panels) for simulations CS (upper panels) and MA (lower 
panels), alongside the evolution of its statistics. The plots show the comoving distance of each satellite from the position of G1/G2 
as a function of redshift, where coloured lines indicate the major (> 30%) and intermediate mass (between 5 and 30%) mergers (see 
Table [2]). The greyscales indicate minor mergers, the tone representing the stellar mass of the progenitor as a fraction of the 
stellar mass of the galaxy at the time of first infall into the corresponding R 200 (see Section |4] of the main text for full details). The 
shaded region indicates 7?2oo of G1/G2, whilst the black horizontal dashed line indicates a radius 30 ckpc, a cut which encompasses 
the baryonic component of the galaxy. The dotted lines follow the positions of the major progenitors after first infall within the 7?2oo 
until subfind can no longer identify them as distinct objects. 


The more violent merger history of G1 is also clear from the 
evolution of /di SC , as G1 exhibits significant morphological evo¬ 
lution between z = 3 and z = 0. Most of the rotational support 
of G1 is acquired between z ~ 2 and z ~ 1.5 in both runs. In 
G1CS, the disc is stable until z ~ 0.3, when it is disrupted co¬ 
inciding with the merger events occurring at z ~ 0.3-0.1. For 
Gl-MA, we find an important decrease in /di SC between z— 1.5 
and z ~ 0.5, and a second phase of disc growth after z ~ 0.4. 
Despite the differences in the evolution of Gl-CS and Gl-MA, 
mainly driven by the different merger times of the satellite galax¬ 
ies (when they enter the inner regions), in both simulations we 
find a galaxy where the rotational support, when present, cannot 
be maintained for long periods of time. 

The more peaceful merger history of the G2 halo is con¬ 
sistent with the rearing of a more disc-dominated galaxy, and 
a more continuous and smooth disc growth. In both simulations, 
G2 evolves into a galaxy with a significant disc, little affected 
by disruptive merger events. An important difference is however 
the typical age and formation time of the discs in the two simu¬ 
lations: the disc of G2-CS is present from very early times, and 
therefore is older on average, while in G2-MA the disc does not 
form until after z = 1 and never reachs /di SC values as large as 
those found for G2-CS. This is, as explained before, the result 
of the stronger feedback assumed in the MA model, which pro¬ 
duces a lower star formation activity at early times due to the 
prompt effects of radiation pressure and kinetic feedback from 
young stars. 


Finally, in the lower panels of Fig. [7] we show the corre¬ 
sponding evolution of the cosine of the angle a between the spe¬ 
cific angular momenta of the gas (j gas ) and of the stars (j starsX 
namely 

COS(cr) = jstars ’ jgas > (5) 

where the hat denotes a normalised vector. In previous work, we 
found that if the angular momentum vectors of the gas and the 
stars are misaligned, the stella r discs become unstable a nd the 
D/T ratios inevitably decrease (IScannapieco et all 120091) . Such 
a misalignment can occur as a consequence of merger events, 
but also if accreting gas brings angular momentum in a different 
direction than that of the pre-existing stellar disc. 

Consistent with our previous analysis, we find that the align¬ 
ment between the gas and stars is not well preserved during the 
evolution of Gl, mainly as a consequence of the major merger 
events of this galaxy, which determine the end points of the 
growth of the discs. The results are very similar for our two 
runs, which is again an indication that results are robust and not 
affected by the particular feedback assumed but rather by the 
merger history. In the case of G2, the gas and the stars are much 
better aligned at all times, in both runs. The most important dif¬ 
ference in this case is at late times, when G2-CS shows a strong 
misalignment, consequence of the mergers that are taking place 
at z < 0.5. As seen before, in G2-MA one of such satellites has 
already been disrupted while the other has not yet reached the 
central object. 
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Table 2. Properties of the major and intermediate mass mergers of G1 and G2 in our two simulations. The Gla-f denotes these 5 
for Gl, and G2a-d the 2 for G2, ordered by their appearance within the 30 ckpc. Z 30 ckpc and Z 200 denote the highest redshift at which 
the most bound subfind particle enters within 30 ckpc and r 2 oo respectively. f^ r ckpc and f^oo denote the subfind stellar mass of the 
satellite as a fraction of that of the main galaxy at the corresponding redshift. f^ 00 is the total mass fraction of the satellite when it 
reaches T^oo- F° r the most massive mergers (Gla, Gld, Gle, G2b) the subfind mass estimate exhibits significant disruption as the 
object enters R 200 , and so for these cases we calculate the ratios 200 Myr earlier to avoid spurious values. In the final column, we 
indicate with 4 sat’ those satellite galaxies that have not merged/disrupted at z = 0. 



ZR200 

ftOt 

JR200 

Run CS 

/r 200 £30 ckpc 

/•star 

J 30 ckpc 

z = 0 

ZR200 

ftot 

JR200 

Run MA 

/r 200 £30 ckpc 

fstar 

J 30 ckpc 

z = 0 

Gla 

2.52 

67% 

51% 

2.33 

21% 


2.52 

121% 

48% 

2.32 

111% 


Gib 

2.35 

5% 

8% 

1.70 

2% 


2.35 

6% 

8% 

2.09 

5% 


Glc 

1.29 

8% 

6% 

0.82 

5% 


1.33 

12% 

10% 

0.81 

14% 


Gld 

0.95 

13% 

14% 

0.32 

5% 

sat 

0.95 

13% 

33% 

0.37 

28% 


Gle 

0.74 

63% 

44% 

0.28 

27% 


0.74 

15% 

80% 

0.08 

21% 

sat 

Gif 

0.72 

64% 

78% 

0.25 

36% 


0.70 

19% 

108% 

0.04 

43% 

sat 

G2a 

1.48 

11% 

11% 

0.81 

5% 


1.48 

13% 

29% 

0.89 

16% 


G2b 

1.44 

9% 

10% 

0.29 

2% 

sat 

1.40 

13% 

33% 

0.44 

9% 


G2c 

1.33 

15% 

28% 

0.60 

9% 


1.33 

18% 

61% 

0.64 

33% 


G2d 

1.14 

7% 

14% 

0.03 

5% 

sat 

1.05 

8% 

37% 

0.02 

9% 

sat 


In summary, by comparing these two histories we see a pic¬ 
ture where the discs are built-up during periods of absence of 
major mergers and alignment between the gaseous and stellar 
discs. In both simulated galaxies the growth of the discs suffers 
termination events when such alignment is disrupted, usually as¬ 
sociated with merging satellites (but can also be produced by 
misaligned gas accretion), either when they pass within 30 ckpc 
or merely have a close approach. 

5. Discussion and conclusions 

We used a simulation of the Local Universe where a Local 
Group-like pair of galaxies with large-scale environment con¬ 
strained to reproduce the environment of the actual LG forms. 
We study how the two most massive galaxies within the simu¬ 
lated LG, candidates for the Milky Way and Andromeda, formed 
and evolv ed. Thejxmstrained nature of the simulation is relevant 
(ICreasev et all 120151) . as it allows the study of the evolution of 
galaxies in an environment similar to our LG, and thus the in¬ 
vestigation of how it affects the properties of the galaxies along 
cosmic time. 

We run two simulations using the same initial conditions, 
but assuming different feedback prescriptions, in order to as¬ 
sess whether our results might in any sense be affected by the 
particular choice of feedback model. It is well known that the 
detailed properties of simulated galaxies depends on the imple¬ 
mentation of the physical processes included, in particular o n 
the way stellar feedback is treated dScannapieco et all 120121) . 
Furthermore, it is still not possible to select the several input pa¬ 
rameters of the simulation from physical principles, and various 
numerical/technical choices can also affect the final outcome of 
simulations. Despite these problems, which will certainly wit¬ 
ness significant progress in the next few years, the predictions 
of simulations in relation to the evolution of properties such as 
morphologies should be still reliable. 

In our simulations, the two main galaxies of the simulated 
LG, Gl (candidate for M31) and G2 (candidate for MW), exhibit 
rather different formation and merger histories although they are, 


at z = 0, both members of the same LG. Gl is much more ac¬ 
tive in terms of the frequency and mass ratios of its mergers 
compared to G2. More importantly, Gl experiences three ma¬ 
jor/intermediate mass mergers at z < 1, contrary to G2 that has 
(since z = 3) only four intermediate mass mergers, and these 
occur before z = 1. The merger rates of the two galaxies are 
then a first indication that the two galaxies will follow differ¬ 
ent evolutionary paths, in particular in relation to the forma¬ 
tion/destruction of stellar discs. 

In fact, we find that Gl, having 6 intermediate/major merg¬ 
ers between z = 3 and z = 0, cannot grow a stellar disc that 
survives until the present time. The various mergers disturb Gl 
significantly, changing the amount and orientation of the angu¬ 
lar momentum of the stars and gas, that in both codes leads to 
a posterior misalignment of the two vectors. When such a mis¬ 
alignment is present, any pre-existing disc gets fully or partially 
destroyed. In contrast, the much more smooth merger history of 
G2, where the angular momentum of the stars and gas are very 
well aligned during most of the time, is consistent with the for¬ 
mation of an extended, rotationally supported disc-like galaxy. 
It is worth noting that mergers, even fairly sizeable ones, do not 
necessarily destroy discs, a nd th e actual damage depends on the 
orbital parameters (see also I Wang et al.ll2012h . 

The different merger/formation histories of Gl and G2 is also 
reflected in their star formation rates and in-situ fractions. Gl has 
a much more violent SFR, and low in-situ fractions, indicating 
a large contribution of stellar material from infalling satellites. 
With such an active merger history, it seems very hard to grow 
a stellar disc with the satellites feeding the disc, instead of de¬ 
stroying it. The more peaceful history of G2 and its much higher 
in-situ fractions at all times suggest that this galaxy will have a 
much higher chance of disc survival. 

We investigated and compared the evolution of the MW and 
M31 candidates in two runs, a first one where only supernova 
feedback is considered (CS) and a second one that additionally 
includes the effects of radiation pressure from massive stars and 
a kinetic component for the supernova energy. Such comparison 
allowed us to assess the robustness of our results and their depen- 
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Fig. 7. The evolution of the disc-to-total measure /di SC of G1 and G2 in our two simulations. The long/short arrows denote the 
entrance of major/intermediate mass mergers in the inner 30 ckpc in the two runs. The lower panels show the evolution of cos (a), 
the angle between the corresponding angular momentum vectors of the gas and stars. 


dence on the assumed feedback model. While it is true that, due 
to the larger feedback effects, the MA model produces galaxies 
with reduced stellar masses compared to CS, both for the pro¬ 
genitor and for the satellites, mergers are as disruptive in the MA 
model as in the CS simulation, where galaxies have in general 
larger stellar masses. In both simulations we find that the stel¬ 
lar disc of G1 suffers several episodes of partial destruction that 
follow merger events, while in G2 a stellar disc grows steadily 
due to its more quiescent evolution. We found (expected) dif¬ 
ferences in the exact merger times of satellites between the two 
runs, not only when they enter R 200 but more importantly when 
they reach the inner galaxy. These can cause some of the satel¬ 
lites to be, at any given time, already disrupted in one run but 
not in the other. More important are the overall differences in the 
star formation rates in terms of the assembly of the stellar com¬ 
ponents: the galaxies in CS are systematically older than in MA, 
which is a direct result of the stronger feedback assumed in the 
latter. 

In summary, by comparing the histories of formation of G1 
and G2, and its relation to the formation/destruction of stellar 
discs along cosmic time, we conclude that the rearing of a disc 
can only be produced during quiet periods where no signifi¬ 
cant mergers and/or close encounters are produced, and when 
the angular momentum vectors of the stars and the gas are well 
aligned. These results are consistent with those found earlier 
in S09, where we simulated the formation of eight galaxies in 
a more isolated environment, giving support to the idea that it 
is the particular formation history of haloes at relatively small 


scales (smaller at least than the LG) that determines their mor¬ 
phological changes along cosmic time, rather than their mem¬ 
bership of a loose group such as the LG. 
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